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ABSTRACT 

Blazars have traditionally been separated into two broad categories based 
upon their optical emission characteristics. Blazars with faint or no emission lines 
are referred to as BL Lacertae objects (BL Lacs), and blazars with prominent, 
broad emission lines are commonly referred to as flat spectrum radio quasars 
(FSRQs). The spectral energy distribution of FSRQs has generally been thought 
of as being more akin to the low-peaked BL Lacs, which exhibit a peak in the 
infrared region of the spectrum, as opposed to high-peaked BL Lacs (HBLs), 
which exhibit a peak in UV/X-ray region of the spectrum. All blazars that are 
currently confirmed as sources of TeV emission fall into the HBL category. Recent 
surveys have found several FSRQs that exhibit spectral properties, particularly 
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the synchrotron peak frequency, similar to HBLs. These objects are potential 
sources of TeV emission according to several models of blazar jet emission and 
the evolution of blazars. Measurements of TeV flux or flux upper limits could 
impact existing theories explaining the links between different blazar types and 
could have a significant impact on our understanding of the nature of objects that 
are capable of TeV emission. In particular, the presence (or absence) of TeV 
emission from FSRQs could confirm (or cast doubt upon) recent evolutionary 
models that expect intermediate objects in a transitionary state between FSRQ 
and BL Lac. The Whipple 10 meter imaging air-Cherenkov gamma-ray telescope 
is well suited for TeV gamma-ray observations. Using the Whipple telescope, we 
have taken data on a small selection of nearby(z<0.1 in most cases), high-peaked 
FSRQs. Although one of the objects, B2 0321+33, showed marginal evidence of 
flaring, no significant emission was detected. The implications of this paucity of 
emission and the derived upper limits are discussed. 

Subject headings: AGN, quasar, VHE, gamma-ray 

1. Introduction 

Blazars, which exhibit the most extreme properties of all active galactic nuclei (AGNs), 
have traditionally been divided into two categories. These categories are separated from one 
another based upon their optical line emission properties. BL Lac objects have either no 
emission lines or weak and narrow emission lines, with the typical definition requiring rest 
frame equivalent widths less than 5 A[Stickel et al. , 1991, Perlman et al. , 1998], while flat 
spectrum radio quasars (FSRQs) exhibit broad optical emission lines. Both classes of objects 
have extreme natures, characterized by a wide range of flux and spectral variability on many 
timescales. Objects from both classes also exhibit high levels of optical polarization. 

The broadband radiation spectrum of blazars consists of two distinct components: a 
synchrotron spectrum component that spans radio to optical-ultraviolet wavelengths (and 
to X-rays for high-frequency peaked objects) and a high-energy component that can extend 
from the X-ray to the very high energy (VHE) gamma-ray (300 GeV < E < 100 TeV) regime. 
Observationally, the spectra appear to have two distinct broad peaks when plotted in a vj v 
representation [Padovani & Giommi , 1995, Fossati et al. , 1998, Ghisellini et al. , 1998]. The 
emission in the first peak of the spectral energy distribution (SED) is generally considered to 
be synchrotron emission from relativistic electrons. The most widely invoked models which 
attempt to explain the higher energy emission of the second SED peak fall into the category 
of leptonic models. These leptonic models posit that the X-ray to gamma-ray emission is 
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produced by inverse Compton scattering of lower energy photons by beamed relativistic 
electrons. The low energy photon fields could arise from synchrotron continuum photons 
within the jet (e.g. Koenigl [1981]), or they could arise from ambient photons from the 
accretion disk which enter the jet directly (e.g. Dermer, Schlickeiser, & Mastichiadis [1992]) 
or after scattering or reprocessing (e.g. Sikora, Begelman, & Rees [1994]). In addition to 
these leptonic models, hadronic models also attempt to explain the second component of 
the SED. These models involve proton-initiated cascades (e.g. Mannheim [1993]) and/or 
proton synchrotron radiation [Muecke and Protheroe , 2001, Aharonian , 2000], as well as 
synchrotron emission from secondary muons and pions [Muecke et al. , 2003]. 

In the blazar sequence proposed by Fossati et al. [1997] and Ghisellini et al. [1998], the 
first SED peak (the synchrotron peak) is at low energies for objects with a high bolometric 
luminosity, and the synchrotron peak is at high energies, sometimes in excess of 100 keV, 
for objects with low bolometric luminosity. In this sequence, FSRQs have their peaks at 
the lowest frequencies, and BL Lac objects have a first peak at higher frequencies. BL Lac 
objects are further divided into two categories: low-peaked BL Lacs (LBLs) exhibit a first 
peak at IR/optical frequencies, and high-peaked BL Lac objects (HBLs) exhibit a first peak 
at UV/X-ray frequencies. 

With the advent of more recent deep X-ray and radio surveys, it is becoming clear 
that there exist objects that do not fall within the confines of this sequence. Some FSRQ 
objects that exhibit intermediate to high frequency spectral peaks, like their BL Lac cousins, 
are being discovered. Padovani, Giommi, & Fiore [1997] created a catalog based on early 
surveys and found that ~17% of the FSRQs in the sample had SEDs similar to HBLs, rather 
than LBLs. Perlman et al. [1998] and Landt et al. [2001] have found that « 30% of 
the FSRQs found in DXRBS have X-ray to radio luminosity ratios characteristic of HBLs, 
while still retaining broad and luminous emission lines characteristic of FSRQs. The recent 
surveys attempt to bridge the gap between previous X-ray and radio selected objects, and 
they indicate the presence of a population of FSRQs with HBL-like SEDs, referred to as 
HFSRQs by Perlman et al. [1998]. The possibility that these HFSRQs could produce 
GeV/TeV gamma-ray emission similar to their HBL cousins is explored in the observations 
described in this paper. 

Although the catalog of confirmed and statistically significant (> 5a) VHE gamma-ray 
sources now includes 8 objects, 6 of which are AGN, the only type of AGN to be decisively 
detected at this time are HBLs [Weekes , 2003] (note: there are 18 VHE objects if one 
includes unconfirmed, yet published, sources). FSRQs with high frequency peaks offer an 
opportunity to expand the VHE gamma-ray catalog and simultaneously constrain the nature 
of HFSRQs, as well as their associated emission mechanisms. By determining the categories 
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of objects that can (and can not) emit VHE gamma rays, the nature of the local medium 
required for such high energy acceleration can be characterized. 

The remainder of this article is organized as follows. Section 2 will identify and describe 
the chosen HFSRQ objects that are candidate GeV/TeV gamma-ray sources. Section 3 will 
describe the Whipple telescope and the observations performed with it. Section 4 will present 
the results of these observations, and Section 5 will discuss the implications of these results. 

2. Candidate Source Properties 

Candidate sources of TeV emission were chosen from lists of FSRQs based upon several 
criteria, including spectral characteristics at radio to X-ray wavelengths, the location of 
the object in the sky, and the redshift of the object. The objects selected for observation 
are shown in Table 1. These candidates were selected from the source lists published by 
Perlman [2000] and Padovani et al. [2002] following the application of several selection 
criteria. Only objects with a high X-ray flux (>10 -12 erg cm" 2 s" 1 in the 0.1-2.4 keV energy 
band) were selected. Objects were also selected by looking at their radio to X-ray spectral 
energy distribution to look for signs of a particularly high frequency of the first SED peak 
when plotted in a vf u representation. The synchrotron peaks for the chosen candidates are 
in the range bounded by ~fewxl0 14 Hz and ~2xl0 16 Hz. All of the sources chosen for this 
study have a high X-ray to radio luminosity ratio, as well as effective spectral indices that are 
more typical of HBLs than classical FSRQs. Another important requirement that was placed 
on the candidate sources was that they be easily observable by the Whipple Observatory, 
which is restricted to observations during clear and moonless nights of Northern Hemisphere 
sources. 

The last restriction that was generally placed on the source candidate list was that they 
be nearby (z<0.1). Sources at higher redshifts begin to have a significant fraction of their 
GeV/TeV emission absorbed due to pair production from the interaction of the high energy 
photon and infrared photons from the diffuse extragalactic background [Gould & Schreder 
, 1967, Vassiliev , 2000]. At this time, the most distant confirmed source of TeV photons is 
at a redshift of z=0.129 (H1426+428), which places it at an optical depth of r^l for 400 
GeV photons [de Jager & Stecker , 2002, Horan et al. , 2002]. All of the candidate sources 
chosen for this study, with the exception of RGB J1629+401, are at redshifts less than z—0.1. 
According to Padovani et al. [2002], RGB J1629+401 has an exceptionally high first peak 
in its SED (~2x 10 16 Hz), so it was included in this survey, despite its redshift of z=0.272. 
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3. Observations 

Observations were performed throughout the 2001/2002 and 2002/2003 seasons at the 
Whipple Observatory, using the Whipple 10-m imaging air Cherenkov telescope (IACT). 

3.1. The Whipple 10-m Telescope and Data Analysis 

The Whipple telescope consists of a 10 meter Davies-Cotton reflector [Davies & Cotton 
, 1957] and a 490 pixel camera [Finley et al. , 2001] comprised of photomultiplier tubes (this 
analysis uses only the inner 379 pixels, which have 0.12° spacing). When a VHE gamma 
ray or cosmic ray hits the top of the atmosphere of the Earth, the resulting cascade of 
particles produces Cherenkov light as it propagates down through the atmosphere. This 
light is collected by the 10-m mirror, and the resulting image that is formed on the camera 
is used to characterize the shower. The parameters used for this characterization are: alpha, 
length, width, distance, size, maxl, and max2 [Hillas , 1998, Punch et al. , 1991]. These 
parameters are described in Table 2, and the geometry of the parameters is shown in Figure 
1. Due to different scattering angles during propagation of electromagnetic showers and 
hadronic showers, as well as the isotropic distribution of cosmic ray shower directions of 
origin, these characteristic parameters can be used to differentiate between gamma-ray and 
hadronic primary particles. By placing cuts on these parameters, the background due to the 
more numerous cosmic rays can be reduced, and thus, the gamma-ray signal-to-noise ratio 
can be increased. These cuts are optimized each season by observing the Crab Nebula, which 
is a steady emitter of VHE gamma rays [Hillas , 1998], and maximizing the signal to noise 
ratio (significance). This procedure, of course, leads to a situation in which the analysis is 
optimized for a Crab-like spectrum. For sources that are significantly softer or harder, one 
has to be resigned to a loss in sensitivity, or one must re-optimize the cuts for a different 
spectral shape. The Crab is also used to calibrate the response of the instrument. The peak 
response of the telescope during the 2001/2002 and 2002/2003 seasons was at ^400 GeV for 
a Crab-like spectrum. 

The Whipple telescope is operated in two standard modes of operation. The first of 
these, referred to as the ON/OFF mode, requires that the 28 minute run in which the tele- 
scope is pointed at the candidate source (ON) be followed by another 28 minute run (OFF) 
in which the telescope is offset in right ascension by 30 minutes such that it images a region 
of sky through the same range of elevation and azimuth as the original source pointing. Since 
this OFF run is taken in nearly the same conditions as the ON run, without the candidate 
source in the field of view, it provides an independent measurement of the background. When 
operating in the other mode, referred to as TRACKING, one observes only at the sky coor- 
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Table 1. Candidate FSRQ Source Properties 



Source 


z 


RA 


Dec 


1 X —ray 


Fradio (mJy) c 


B2 0321+33 


0.062 


03 24 41.1 


+34 10 46 


6.6 a 


368 


PG 2209+184 


0.070 


22 11 53.8 


+ 18 41 52 


8.4 a 


116 


WGA J0838+2453 


0.028 


08 38 11.1 


+24 53 45 


6.5 a 


32 d 


RGB J1413+436 


0.090 


14 13 43.7 


+43 39 45 


4.5 b 


50 


RGB J1629+4008 


0.272 


16 29 01.3 


+40 08 00 


9.0 b 


20 



a 0.3-3.5 keV flux (10 12 erg cm 2 s x ) corrected for galactic absorption taken 
from Perlman [2000] 

b 0. 1-2.4 keV flux (10~ 12 erg cm -2 s^ 1 ) uncorrected for galactic absorption 
taken from Perlman [2000] and Laurent-Muehleisen et al. [1998] 

c 4.85 GHz radio flux taken from Gregory & Condon [1991] 

c 4.85 GHz radio flux taken from Condon, Anderson, & Broderick [1995] 



Table 2. Image Parameters Used for Characterization of Showers 



Parameter 


Description 


length 


length (major axis) of shower image 


width 


width (minor axis) of shower image 


distance 


distance from shower image centroid to camera center 


size 


sum of signals from all pixels in shower image 


maxl 


largest signal recorded in any image pixel 


max2 


second largest signal recorded in any image pixel 


alpha 


angle of shower major axis relative to line from camera center to image centroid 
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Fig. 1. — Image characterization of air showers imaged at the focal plane of the telescope. 
This cartoon shower projection has an elliptical shape that is typical of gamma ray and 
cosmic ray induced showers, and it has a large alpha angle, which implies that the field of 
view is not centered on the source. 
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dinates of the source. This mode utilizes the large data set of OFF source observations taken 
throughout the season to estimate the background expected at small alpha angles (where the 
source counts are expected for point-like sources), as a function of the background counts at 
large alpha angles (where no source counts are expected). This procedure takes advantage 
of the fact that gamma-ray events from the source will preferentially produce an image with 
a small alpha angle, while the more numerous background cosmic rays will be randomly 
oriented in the camera field of view. While this method allows more on-source data to be 
taken, it requires an accurate calculation of the tracking ratio, which is the ratio of events 
with alpha angle between 20° and 65° to events with alpha angle between 0° and 15° in 
the absence of a gamma-ray source. Since the analysis results are highly dependent on this 
tracking ratio, it must be calculated independently each season using OFF region data that 
have been taken when the telescope is operating in a state that is as close as possible to the 
state of operation when the TRACKING data is taken. This ratio has been calculated and 
applied to this analysis. For the 2001-2002 data, the tracking ratio is 0.312 ± 0.003; and for 
the 2002-2003 data, the tracking ratio is 0.3067 ± 0.0035. 

Once the tracking ratio has been calculated using OFF source observations, the back- 
ground for the ON source observations in the 0° to 15° alpha region, where all of the gamma- 
ray signal is expected, can be calculated. This is done by using the tracking ratio to scale 
the background events in the 20° to 65° alpha region to obtain the expected background 
counts in the 0° to 15° alpha region. After applying a simple propagation of error formula 
to this procedure, one obtains the following for the significance. 

_ Nqo_ 15 o — r(N20°-65°) qx 

vW°-i 5 °) +r 2 A 20 o_ 65 o + (Ar) 2 (jV 20 o_ 65 oy' 

where r ± Ar is the tracking ratio and its associated statistical error. An important char- 
acteristic of this method for computing significance is that it incorporates the error on the 
tracking ratio, which is a quantity that is not accommodated by the method of Li & Ma 
[1983]. It was found by Li & Ma [1983] that this type of significance calculation tends to 
be more conservative than their reported method of calculating significance. It should be 
noted that the error used to calculate this significance value is the statistical error; it does 
not include possible systematic errors. 

A more detailed description of the standard analysis and the use of various cuts has 
been published previously in several articles including Mohanty et al. [1998], Reynolds et 
al. [1993], and Catanese et al. [1998]. 
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3.2. Source Observations 

Each of the candidate FSRQs was observed by the Whipple telescope, as shown in Table 
3. Two of the candidates were observed during the 2001/2002 season (B2 0321+33 and PG 
2209+184), and three additional candidates were observed during the 2002/2003 season 
(WGA J0838+2453, RGB J1413+436, and RGB J1629+4008). One source, B2 0321+33, 
was observed during both seasons since there was a marginal hint of flaring during the 
first season that merited a follow-up set of observations. During both seasons, the image 
parameter cuts were optimized independently, but the optimized cuts were identical for the 
two seasons, as expected for a stable camera configuration. Tracking ratios were calculated 
from a large number of OFF source runs taken during the relevant time frames. All of the 
observations reported here were taken at zenith angles less than 35°. 

4. Results 

Significance values calculated using the method described above are shown in Table 3. 
None of the sources exhibited significant steady emission of gamma rays over the integrated 
time frame of the observations. Also shown in Table 3 are the 95% upper limit values for 
steady emission, calculated using the method of Helene [1983]. The fifth and sixth columns 
of Table 3 are the upper limits for flux detected at Earth, and the last column shows the 
upper limits for the flux at the source. This last column has been calculated using the optical 
depth of the extragalactic background for 400 GeV photons, as given by the baseline model 
of de Jager & Stecker [2002]. 

Although none of the FSRQs showed evidence for steady emission during the extended 
observations, one of them did exhibit marginal evidence of flaring during the 2001/2002 
observing season. On 22 Oct. 2001 (MJD 52204) two 28 minute observation runs were taken 
on B2 0321+33. The average gamma-ray rate of these two runs was 0.46 ± 0.14 times that 
of the Crab, which corresponds to a significance of 3.3cr. The peak rate was 0.62 ± 0.19 
Crabs. The light curve of data taken during that season is shown in Figure 2. If the trials 
factor is derived using the 16 observation nights of FSRQs during this season (combined for 
all sources), then the post-trials significance of this rate increase is 2.5a. While this rate 
increase was not significant enough to claim a detection of a variable source, it was enough 
to merit the further observations that were taken of this candidate in the 2002/2003 season. 
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Fig. 2. — Light curve of B2 0321+33 using all data taken in the 2001/2002 observing season. 
Error bars are la statistical errors. 
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5. Discussion and Conclusions 

No significant emission has been detected from any of the candidate sources in this 
initial survey. There was marginal evidence of a rate increase observed in the B2 0321+33 
light curve, but the statistical significance of this increase is 2.5a (post-trials significance), 
which could be accounted for by a statistical fluctuation. None of the other objects showed 
significant steady emission or any interesting features in their light curves. 

The upper limits presented here are the first VHE gamma-ray upper limits published 
for HFSRQs; however, they do not severely constrain the emission models due to the large 
amount of uncertainty in the models (and, of course, due to limited sensitivity of the tele- 
scope). Potential variability of the sources also leads to uncertainty in the expected flux. 
Typically, FSRQs and BL Lacs are highly variable sources. Any VHE gamma-ray flux de- 
tected from these objects would most likely be from the sum of many flares or from a large 
flare that happened to be caught during the time frame of the observation. Applicable mod- 
els of FSRQs have been made by Padovani et al. [2002]. In this paper, several high-peaked 
and intermediate HFSRQs were modelled using a one-zone synchrotron inverse Compton 
model that was developed by Ghisellini, Celotti, & Costamante [2002]. This model includes 
synchrotron self-Compton, as well as seed photon contributions from the disk and the broad 
line region. The inclusion of the latter components may be important for these sources, 
which are known to have broad line emission. While Perlman [2000] predicted that B2 
0321+33, PG 2209+184, and RGB 1413+436 would be good TeV emission candidates, no 
modelling of the TeV flux was reported. 

As a rough estimate of the expected flux, we have compared the SED of B2 0321+33 to 
those modelled by Padovani et al. [2002], and we find that it is similar to WGA J0546.6- 
6415, which has an estimated intrinsic very high energy gamma-ray flux of ~ 2 x 10~ 12 erg 
cm~ 2 s _1 above 400 GeV. This is below the value of our reported upper limit. With the 
exception of RGB J1629+4008, the very high energy flux of the other objects is not as easy 
to predict within the context of published models. More modelling is necessary. 

Padovani et al. [2002] directly modelled the SED of RGB J1629+4008, and they derived 
a synchrotron peak at 1.7 x 10 16 Hz and a very high energy gamma-ray flux of ~ 2 x 10~ 14 
erg cm -2 s" 1 above 400 GeV. No mention was made of absorption of gamma rays by the 
diffuse IR background. This will produce a significant modification since the source is at 
a large redshift of z=0.272. Using the model of de Jager & Stecker [2002], this redshift 
corresponds to an optical depth of sal. 8 for 400 GeV gamma rays. The intrinsic source flux 
upper limit derived here is 0.29 x 10~ 12 erg cur 2 s _1 above 400 GeV, which is certainly 
higher than the predicted source flux so no strong constraints can be placed on the steady 
emission models for this object. The particularly high synchrotron peak and the potential 
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for increased flux during flaring prompted these observations, but no variability was detected 
during the observations. 

Cavaliere & D'Elia [2002] and Bottcher & Dermer [2002] posit that the blazar sequence 
(i.e. the relationship of blazar type and synchrotron peak location to the bolometric lumi- 
nosity) is fundamentally dependent on the accretion rate, m. In this scenario, objects with 
large values of rii are FSRQs and have large Lft / (more than ~10 46 erg s _1 ); while objects 
with small values of hi are BL Lacs and have small (~10 44 erg s _1 ). Both of these 
models also predict an evolutionary scenario that progresses from FSRQs to BL Lacs, which 
leads to intermediate FSRQ objects that can have optical lines and large L^, along with a 
synchrotron peak more like that of BL Lacs. These objects should be excellent candidates 
for VHE gamma-ray emission. The fact that the observations presented here yield little or 
no VHE emission does not support this scenario, but these observations certainly do not 
rule out the proposed models for two reasons. It is possible that these intermediate objects 
have optical depths that are large enough to prevent the escape of VHE emission. It is also 
possible that the handful of objects in this initial survey were not in a high flaring state 
at the time of observation. The inherent variability of blazars makes the detection of both 
BL Lac and FSRQ VHE emission difficult, as evidenced by the paucity of detected sources 
relative to the number of observed BL Lac candidates in past campaigns [de la Calle Perez 
et al. , 2003, Horan et al. , 2004]. More observations of these and other objects will be 
necessary. 

Although no detection of VHE gamma-ray emission has been reported from this initial 
survey, dim or variable emission from these candidates, or emission from other similar candi- 
dates, is not ruled out. This work should be continued with more sensitive instruments, such 
as the next generation of IACTs (VERITAS, HESS, MAGIC, and CANGAROO-III) that are 
beginning to come on line at the time of this writing [Krennrich et al. , 2004, Hinton et al. 
, 2004, Lorenz et al. , 2004, Kubo et al. , 2004]. More modelling of these objects is required 
in order to predict the high energy flux more accurately, and therefore, the potential of a 
given HFSRQ to be a gamma-ray source. In addition, the continued progress of deeper X-ray 
to gamma-ray surveys (Swift, EXIST, INTEGRAL, GLAST) and X-ray imaging (Chandra, 
XMM) may reveal additional TeV FSRQ candidates with higher frequency peaks in the SED 
and/or increased high energy flux potential. 

We acknowledge the technical assistance of E. Roache and J. Melnick. This research is 
supported by grants from the U.S. Department of Energy, Enterprise Ireland, and PPARC 
in the UK. 
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Table 3. Whipple Observations of Candidate FSRQ Sources 



Source 


ON2001-2002 
(min) a 


ON2002-2003 
(min) a 




upper limit c 
(Crab) 


upper limit d 
(10- 11 erg cm- 2 s" 1 ) 


upper limit 
(10- 11 erg cm" 2 s- 1 ) 


B2 0321+33 


320.3 


331.5 


1.56 


0.10 


0.52 


0.68 


PG 2209+184 


217.6 





0.15 


0.13 


0.71 


0.97 


WGA J0838+2453 





781.6 


0.15 


0.05 


0.29 


0.32 


RGB J1413+436 





535.4 


0.75 


0.06 


0.35 


0.53 


RGB J1629+4008 





415.1 


1.29 


0.09 


0.47 


2.9 



a ON refers to minutes spent on source during stable/clear weather conditions 
b significance calculated as described in text 

c 95% confidence level upper limits calculated using method of Helene [1983], expressed as a fraction of the Crab flux (Rcrab = 
2.84 ± 0.16 in 2001-2002 season and R Cra b = 2.55 ± 0.13 in 2002-2003 season) 

d 95% confidence level upper limits for emission above Ri400 GeV, assuming a Crab-like spectrum as in Hillas [1998] ; calculated 
using method of Helene [1983] 

e 95% confidence level upper limits as described above, including the effect of absorption from the cxtragalactic background, using 
the baseline model of de Jager & Stecker [2002] 



